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Abstract 

The analysis of spatial extremes requires the joint modeling of a spatial 
process at a large number of stations and max-stable processes have been 
developed as a class of stochastic processes suitable for studying spatial 
extremes. Spatial dependence structure in the extreme value analysis 
can be measured by max-stable processes. However, there have been few 
works on the threshold approach of max-stable processes. 

We propose a threshold version of max-stable process estimation and 
we apply the pairwise composite likelihood method by [H] to estimate 
spatial dependence parameters. It is of interest to establish limit behav- 
ior of the estimates based on the settings of increasing domain asymptotics 
with stochastic sampling design. Two different types of asymptotic nor- 
mality are drawn under the second-order regular variation condition for 
the distribution satisfying the domain of attraction. The theoretical prop- 
erty of dependence parameter estimators in limiting sense is implemented 
by simulation and a choice of optimal threshold is discussed in this paper. 

1 Introduction 



Extreme value theory and its application are dealing with related methodologies 
to understand phenomena of rare events such as flooding, high temperatures 
and precipitations in environmental data. The behavior of rare events requires 
understanding of the tail distribution. 

Extreme value theory has been studied for the univariate case in which 
extremes are observed as a single variable, during a few decades since [TT] and 
[TB] . p5] and \T provides statistical methods in the analysis of extremes, and [3] 
is a very useful reference with the introduction of modeling and applications of 
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extreme values. Multivariate extreme value theory has been developed to build 
the modeling of joint extremal behavior. |21j reviewed relevant theories in the 
view of probability and measure theory for multivariate extremes. 

In a spatial context, a single quantity (e.g., sea level) is measured at multiple 
locations and the observed data are spatial variables which are distributed across 
the earth's surface. Therefore one ultimately requires the modeling of spatial 
extremes, and a spatial dependence among the different locations is of interest. 
[B] introduces several references with issues of spatial extremes. 

It is natural to consider a stochastic process when the sample maxima are 
observed at each site of a spatial process. Max-stable processes have been 
developed as an infinite dimensional generalization of multivariate extremes. 
The first general characterization of max-stable processes was by [7], and [27] 
has constructed a special case of max-stable processes which provides the useful 
interpretation of extreme rainfall models. Statistical techniques based on the 
Smith's max-stable model have been developed by [3] and [S] and the well-known 
classes of max-stable processes are discussed further by [23] and [13]. However 
the modeling of max-stable processes did not give a straightforward usage due 
to the complexity and unavailability of the full likelihood for the max-stable 
model, and [19) developed the maximum composite likelihood approach to fit 
max-stable processes. 

Though max-stable processes for blocked maxima approach are on the ex- 
ploratory stage, the research on max-stable processes with exceedances over 
threshold has hardly been considered. In this paper, we are concerned with 
the development of a threshold approach using max-stable processes in spa- 
tial extremes. We review the background of extreme value theory, max-stable 
processes and spatial dependence measure in Section [2] In Section [3| we intro- 
duce our methodology to model exceedances over threshold using max-stable 
processes and describe its theoretical framework. Section |4| develops asymptotic 
properties of spatial dependence parameter estimates, which are illustrated with 
a simulation study. 

2 Modeling of spatial extremes 
2.1 Extreme value theory 

Let Xi, • • • , Xn be i.i.d. random variables with the same probability distribution 
F and let Af„ = max(Xi, • • • be the maximum. If il/„ converges under 

renormalization to some nondegenerate limit, then the limit must be a member 
of the parametric family, i.e. there exist suitable normalizing constants a„ > 0, 
bn and the distribution G such that 

p{ ^" ~ ^" <x\^ F'^ianx + b„) G{x), asn^oo (1) 
L a„ J 

where G is a nondegenerate distribution function. The distribution function 
G which is possible limit laws for maxima of i.i.d. sequences has one of three 
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Extreme Value Distributions (EVD). 

The three types of EVD can be represented as G combining into a single 
parametric family distribution, which is called the Generalized Extreme Value 
(GEV) distribution: 

G{x; /X, V, = exp I - (^1 + ^^^^ 

where y+ = max(0, y), ^ is a location parameter, cr > is a scale parameter and 
^ is a shape parameter which determines the tail behavior. The Generalized 
Extreme Value distribution G has a max-stable property: ii Xi, - ■ ■ ,Xn are 
i.i.d. from G, then max(Xi, • • • ,Xn) also has the same distribution, i.e. 

G'^{x) — G{Apjx + Bpj) for existing constants > 0, B^. 

|16| showed the relationship between extreme value distributions and max-stable 
distributions that any extreme value distribution is max-stable and vice versa. 

The form of the limiting distribution is invariant under monotonic trans- 
formation. Therefore, without loss of generality we can transform the GEV 
distribution into a specific standard form, called unit Frechet distribution, 

p| (l + ^Mll—l^) z I = P{Z <z)= exp(-l/z), z > 0, 



and note that the unit Frechet form is a distribution which has the max-stable 
property. 

Multivariate extreme value theory is concerned with the joint distribution 
of extremes of two or more random variables. If G is a multivariate EVD, the 
marginal distribution must be represented by the GEV distribution and each 
marginal GEV distribution can be transformed into unit Frechet margin, which 
has the max-stable property. 

The finite-dimensional framework of multivariate extreme distribution is ex- 
tended to an infinite-dimensional generalization with spatial processes. The 
infinite-dimensional extremes has quite analogous extension to the theory of 
max-stable random vector. Let 5 be a study region and denote s as a location 
in the study region. If there exist normalizing sequences a„(s) and bn{s) for all 
s G S such that the sequence of stochastic processes 

nrax ^M-^^Y{s) (2) 

where Y(s) is non-degenerate for all s, then the limit process Y{s) is a max- 
stable process. A finite sample {Y{si), • • • , Y{sd)} can be concerned as a real- 
ization of a spatial process Y{s) for more realistic setting. 

2.2 Max-stable process and composite likelihood 

Suppose X{s),s & S is a stochastic process, where 5 C R'^ is an arbitrary index 
set. We can interpret X{-) as a spatial process and an appropriate generalization 
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of multivariate extremes can be made in terms of spacial processes as following: 
for each n > 1, there exist continuous functions a„(s) positive and 6„(s) real, 
for s S 5 such that 

p^nf Xis,)-b is,) ^ ^^^^^^^ ^ 1^ . . . j _^ G,,,.. ,sA^i^l): ■ ■ ■ ,^{sk)). 

(, O-n [Sj ) ) 

_(3) 

Then Gs-i,---,sk is £^ multivariate extreme value distribution and the limiting 
process is max-stable if ([3]) holds for all possible subsets si, • • • , s/f G 5. Note 
that this is equivalent to the expression in equation 

Wc are interested in modeling and estimation using max-stable processes for 
extremes observed at each site of a spatial process. A general representation 
of max-stable processes was first given by The conceptual idea of max- 
stable processes can be constructed by two components: a stochastic process 
W{s) and a Poisson process 11 with intensity dC/C^ on (0, oo). If {W^i(s)}igN is 
independent copies of W{s) with E\W{s)] = 1 for all s and Ci G H, i > 1, is 
points of the Poisson process, then 

F(s) = max CiWi(s), s€S 

i>l 

is a max-stable process with unit Prechet margins. The construction of different 
max-stable processes can be differentiated from different choices of the W{s) 
process and the well-known classes of max-stable processes are discussed by 
[27], [24] and [M]. 

Bivariate joint distribution is derived for each max-stable process model. 
The Smith model has the exact form of bivariate distribution 

P{Yisi)<yi,Yis2) < 

= exp i - — $ ^ + - log ^ I - — $ I ^ + - log ^ I [. (4) 




where a = \J [s\ — S2)^Tj~^[s\ — s-i) characterizing spatial dependence with co- 
variance matrix S and <& is the standard normal cumulative distribution func- 
tion. Bivariate marginal distribution of Schlather model is given by 









(-51 







e-P<! - + + 2(P(M + 1) ^^I'll^y I ) (5) 

and the correlation p(K) represents spatial dependence where h is the Euclidean 
distance, ||si — S2II, between two stations. Max-stable model of 14 is called 
the Brown-Resnick process and the closed form of the bivariate distributions 
associated to the variogram 7 is given by 
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Pairwise composite likelihood We are interested with the analysis of spatial 
extremes at a large number of stations and the standard methods of estimation, 
such as MLE and Bayes methods, require a full likelihood. However the full like- 
lihood for the max-stable processes may not be available due to the complexity 
of its analytic form. With the lack of an explicit form of the joint distribu- 
tion, [IH] developed a pairwise composite likelihood approach to fit max-stable 
processes, based on a composite likelihood method by |17) . 

Assume M i.i.d. replications of a stochastic process with bivariate densities 
fiVii yjj4')j 1 1^ hj ^ K^ina spatial region with K locations. Then the pairwise 
composite log-likelihood is defined by 

M K-1 K 

^•p(V';Y)=^^ ^ w^j\ogf{yrm,ymj;ip), (7) 

m— 1 i—1 

where (i,j) is a pair of stations and is nonnegative weight functions. One 
may set the weight as an indicator function, i.e., = 1 if || si — S2 || < (5, and 
otherwise. The maximum pairwise composite likelihood estimator (MCLE), ■(/;, is 
chosen to maximize [19] stated the asymptotic properties of MCLE based on 
the joint estimation, which maximizes the pairwise composite likelihood instead 
of the full likelihood. 

2.3 Dependence of spatial extremes: Extremal coefficient 

In the analysis of spatial extremes, one can be interested with measuring spatial 
dependence among locations and a metric characterizing the tail dependence 
is extremal coefficient. Suppose a d-dimensional random variable X has the 
common marginal distributions F(x). The extremal coefficient 9d can be defined 
by the relation 

Pr{max(Xi, ■■■ ,Xd) <x) ^ F^\x). 

Assuming the standard form of unit Frechet distribution on each margin, we can 
characterize the dependence among the components of marginal distribution 
independently. Let Z be c?-dimensional maxima with unit Frechet margins and 
whose multivariate extreme value distribution is expressed as 

Pr{Zi <zi,---,Zd<Zd} = cxp{-y(zi, • • • , Zd)}, (8) 

where the exponent measure y is a homogeneous function of order —1. Due to 
the homogeneity of V, the extremal dependence can be measured by V which 
implies complete dependence if V{zi, • • • , Zd) = max • • • , ^) and complete 

independence if V{zi, ■ ■ ■ ^ Zd) — + ■ ■ ■ + j^- 

The relationship between the extremal coefficient 6d and the exponent mea- 
sure V is drawn from 9^ — V{1t--- ,1), and ([5]) is expressed in terms of extremal 
coefficient 

Pr{Zi<z,--- ,Zd<z} = exp(-— 1, (9) 
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where 1 < Od < d with the lower and upper bounds corresponding to complete 
dependence and complete independence, respectively. 

We consider a pairwise extremal coefficient as a special case of © in the 
spatial domain. Let Y{s) be a spatial process with unit Frechet margin for all 
s € S and then extremal dependence between different sites s and s' is obtained 

by, 

Pr{Yis)<y,Y{s')<y} = exp 

A naive estimator of the pairwise extremal coefficient is proposed by |27| . [2 5) 
investigated theoretical properties of the extremal coefficients and proposed self- 
consistent estimators of 9 (i.e. estimators that satisfy the properties of extremal 
coefficients) for the multivariate and spatial case. 



0{s - s') 



3 Modeling for exceedances over threshold 

Consider the distribution of all observations X over a high threshold u and let 
Y = X - u> 0, then 

F„(y)^Mr<.|r>o}= ^(V'l7f"^ - 

I — r (u) 

As u xq = sup{a; : F{x) < 1}, we can find a limit H called Generalized 
Pareto Distribution (GPD) 

Fuiy)^Hiy;a^,0^1- (l+C^^ (10) 

[26] described the bias versus variance tradeoff in the choice of threshold u of 
univariate case. If the threshold u increases, the variance of estimators will be 
high due to small N (number of exceedances) while the estimates are biased 
due to the poor approximation of Fu{-) by H{-) if u is too small. Thus limit 
theorems on the threshold approach in the literature are presented as — ^ oo 
and u = UN ^ xq simultaneously. 

PU] established the rigorous connection between the classical extreme value 
theory and the generalized Pareto distribution and proved that the limit of the 
form (jlOp exists if and only if there exist normalizing constants and the limiting 
form of H such that the classical extreme value limit ([T]) holds. Thus the limit 
result for exceedances over thresholds is equivalent to the limit distribution for 
maxima in this sense. 

As in the univariate case, the threshold method has been developed in 
bivariate case as well. Let {xo,yo) denote the upper endpoint of F, where 
{xo,yo) = sup{(x,?/) : F{x,y) < 1}, and define the conditional distribution of 
{X — u.Y — v) given X > u or Y > v, 

F (xv)- Fi^ + ^^^ + y)-Pi^'^) 
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Then the conditional distribution of bivariate exceedances converges to H where 
H is a. multivariate generalized Pareto distribution by [23]. 

In this section we develop an alternative methodology for threshold ex- 
ceedances using max-stable processes with unit Frechet margins. We suggest the 
modeling of the bivariate threshold exceedances by assuming that the asymp- 
totic distribution holds exactly above a threshold and it leads to a simplified 
dependence structure for max-stable processes as we characterize the depen- 
dence among the components of bivariate marginal distribution in Q , (O and 



The likelihood representation for this threshold method is developed to fit the 
model and this has a similar idea by [22] which establishes a joint distribution for 
Markov chains where the bivariate distributions were assumed to be of bivariate 
extreme value distribution form above a threshold. The censored threshold- 
based likelihood approach is also available for the modeling of spatio-temporal 
extremes in [T5] . 

Threshold methodology Suppose we have annual maxima {i^cs} at site s = 
1, • • • ,D in year t* = 1, - ■ ■ ,T* . We assume the vectors {^t's} are independent 
for different t* with joint densities given by a max-stable process, i.e., an explicit 
expression for its bivariate joint distribution is known and the marginal distri- 
butions are unit Frechet for each t* and s. Then the joint bivariate distribution 
of the annual maxima, Fam is written by 



where 9 is the dependence parameter which can be estimated by the max-stable 
model. Now suppose that the daily data are {Xts, t — I,-- - ,T,s — I,-- - , D} 
and the joint bivariate distribution function is FoAixts,Xts'',0). Assume that 
the daily data Xts form i.i.d. random processes and the annual maxima are 
Yfs. Then the relationship between their bivariate distributions is 



where M is the number of days in a year. We can have a closed form for Fam 
from the max-stable theory and also get an expression for Foa from the above 
representation. 

In practice, we would expect to apply some notion of thresholding. Suppose 
we fix the threshold u and we assume that the same threshold for all locations 
for convenience. Then we observe exceedances {Xts} such that Xts > u. Let 
5s = I{Xts > u) where / is the indicator function. We can obtain the following 
joint distribution of {Ss-, Xts-, ^s', Xts') from four possible regions by including or 
excluding the interval over threshold u, 



FAM{yt*siVt*s']0) ^Pr{Yt*s < yt's^Yt's' < yt's';&}, 



FDA{xts,Xts';0) = Pr{Xts < Xts, Xts' < xts';0} 

^FAM{xts,Xts';ef'''' 



(12) 



Pr{5s = 0, 5s. - 0} 
Pr{5s - l,5s' =Q,Xts < Xts} 
Pr{5s^Q,5s. = I, Xts' < Xts'} 
Pr{5s = 1,5s' = 1} 



Fda{u,u) 

FDA{xts,u) 
FDA{u,Xts') 
FDA{xts,Xts')- 
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We extend the threshold version of max-stable processes and apply the max- 
imum composite likelihood method on it. The likelihood contribution of the pair 
{xts , Xts' ) derived from the joint bivariate density can be obtained by 



L{Xt,,Xt,r,e,rj) = 



Fda{u, u) if Xts < u, Xts' < u, 

g§^FDAixts, u) if Xts > U, Xts' < u, 

■g^FDA{u,Xts') if Xts < U,Xts' > u, 

FDA{xts,Xts') if Xts > U.Xts' > U. 



where 6 is the dependence parameter vector and r/ is a vector of marginal 
GEV parameter. Combining the above likelihood representation with a pairwise 
likelihood, we assume T i.i.d. replications of a stochastic process with bivariate 
densities of the unit Frechet margins L{Xts, Xts' ',d,r]),l < s,s' < D. Then the 
pairwise composite log-likelihood for a thresholded process is 

T D-1 D T 

l{d,v) = Y.Yl E Wss'\ogLiXts,Xts';9,q) = J2ltie,v) (13) 

t=l s=l s' = s+l t=l 

where lt{0,T]) = Y^fji Y.f'=s+i'^'^'>' ^'^sH^ts, Xts'] 6,110), (s, s') is a pair of 
different stations and T is a number of observations. In practice, the marginal 
parameter t] will be estimated but we let r] be the true value r/g to simplify 
theoretical justification. Thus we fix the marginal GEV parameters rj = ijo 
and estimate the dependence parameter 9. A dependence parameter 6 can be 
estimated by maximizing the pairwise composite likelihood function (jl3p with 
the known value r/Q. 

Suppose X^*) — (Xts, Xts') and denote the composite score functions by 
pairwise log-likelihood derivatives as 

^ ' ' de 

t=i 

Then the estimating equations 

;XW,...,X(^))=D,Z(0,r,o;XW,...,X(^)) = O. 



The parameter estimator is a root to solve above estimating equations and 
we now start to describe the theoretical framework with more strict conditions 
to obtain asymptotic properties of the estimator. 



3.1 Second-order regular variation 

To obtain a limiting distribution of F^^v we assume a strict form of condition, so 
called the second- order regular variation condition, for the distribution satisfying 
the domain of attraction. The ideas of second-order regular variation have been 
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applied to the statistics of extremes. Asymptotic properties of estimators in 
univariate extreme value theory have been investigated with the second-order 
regular variation (see [55], [5], and [10]), and the second-order regular variation 
condition was studied for bivariate extremes by [S]. 

Definition 1. A function f{x) is regular varying with index ti if for some 

Ti e R, 

t^oo f{t) 

The function f{x) is second- order regular varying with the first order ti and the 
second order T2 if there exists a function q{t) ^ as t — > cx) such that 

/(to) 



lim , , = x''\ x>0. 

t^oo qi^t) 

Just as in the univariate case, the representation of bivariate regular variation 
exists. 

Definition 2. A function f{x,y) : R"^ — )• is regular varying of index t if 

f{tx,ty) 

t-+oo j[t,t) 

where r(Xx,Xy) = X^r^x^y) for some A > 0. 

See [22] for the related discussion of multivariate regular variation. 

Suppose that {Xi,Yi), i = 1, • • • , n, is a sequence of i.i.d. random vectors 
and F be the common distribution of {Xi,Yi) with marginal distributions Fi 
and ^2- A distribution function F is said to be in the domain of attraction of 
a distribution function G, shortly F G F){G), if 

lim F"{anX + b^^Cny + d„) = G{x, y), a„, c„ > and 6„, d„ e R (14) 

n— ^00 

for all X and y. The two marginals of G{x, 00) and G{oo, y) are one-dimensional 
extreme value distributions satisfying 

lim Fi"(a„a; + 6„) = exp{-il + (ixy^^^'}, 

n— ^00 

lim F2"(c„y + d„) - cxp{-(l + 6j/)"'/«^} 

n— f 00 

where ^1 and ^2 sue real parameters. 

Let (a;o,2/o) denote the upper endpoint of F{x,y) and the conditional dis- 
tribution of {X ~ u^Y ~ v) given X > m or F > w is defined as in pT]) . The 
equation (1141) by taking logarithms can be expressed as 

lim t{l^ F{atx + ht,cty + dt)] = - log G(.t, y) =: $(x,y) (15) 
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and it is checked easily that (ITSI) implies that 

p , , r f, t{l-F{atx + bt,Cty + dt)} 
Im Fb,^dAatX,Cty) = hm 1 -y- -. 

- log G(a:, 7/) 

where is a bivariate generalized Pareto distribution. It has been illustrated 
that H is 8L good approximation of -Ffcj.df in the sense that 

lim sup \Fbt,dt{atX,cty) - H{x,y)\ =0, 

t~¥oo Q^(^atX,Cty)<{xo-bt,yo-d.t) 

if and only if F is in the maximum domain of attraction of the corresponding 
extreme value distribution G. 

Suppose that the second-order regular variation condition in ^ holds: there 
exists a positive or negative function a with lim(_>.oo ct{t) =0 and a function Q 
not a multiple of $ such that 

t{l - Fjatx + bt,cty + dt)} ~ <^{x,y) 

hm -i^ — i = Q{x,y) (16) 

t-^co a[t) 

locally uniformly for {x,y) G (0, cx)] x (0, cx)]. Define Ui as the inverse function 
of 1/(1 — Fi), i = 1,2 and it is known that for x,y > 0, 

UUt'x) -UUt') x^i-l 
lim 7777 = , 

t'-i-oo o,[t ) 7l 

,. U2{t'y) - U2{t') y^- - 1 
lim — = . 

t'-i-oo c(r) 72 

For t 7^ t', define at, &(, Ct and dt such that a(t') = at,Ui{t') = bt,c{t') = ct, 
and U2{t') = dt respectively. Let 

Ui{tx)-bt 

xt 



Vt 



at 

U2{ty)-dt 



Ct 

and we could rewrite the form as 

hm t{l - F(C/i(te), C/2(i2/))} = -logG , ^ <i>o(x,2/). 

It follows the similar form of the second-order condition (|T6)) . 

l--F('Ji( i-f"(,,) ),e/2( i_p^(^^) )) _ 

hm ^£(^iA) ^ ^J^^Q(-Jl^y:i^\, (17) 



,1-F((,t,<it)^ 



71 72 



We can rewrite the condition ([T7| and the following second order condition holds 
for Fb,.df 
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Condition. There exists a positive or negative function A{-) such that 



Fbt,dt{atx,cty) = H{x,y) + A{t)^{x,y) + Rt{x,y), for all t and x,y > (18) 
(i) * = 0, A{t) = o(l) and Rt{x,y) = o{A{t)) as t ^ oo, or 

(a) ^ is continuous and not a multiple of H, A{t) = o(l) and Rt{x,y) = 
o{A{t)) as t^ oo. 

The second order regular variation condition implements tlic domain of at- 
traction condition as a special asymptotic expansion of the conditional distri- 
bution Fb^^dt near infinity. The asymptotic behavior of tail distribution turns 
out to depend on how the regular variation condition behaves. 

In order to obtain asymptotic properties for 6, we need to understand the 
behavior of ED(0o) given the second-order regular variation, where D is the score 
functions of pairwise composite likelihood. The following defines the statement 
on how integrals of the score functions behave corresponding to the second-order 
condition. 

Proposition 1. Let gt{x,y) be any measurable function. Suppose Fb^^dt satisfies 
the condition with (i) or (ii) with function A. Define fbt,dt = dxdv* ' ^(^^v) ~ 



'^andi,{x,y) = qSlM,if 



9t{x,y) 



fbt,dt{atx,cty) - h{x,y) 
A{t) 



il){x,y) 



<K{x,y) 



(19) 



which K{x,y) is integrable, then in case of (i) 

/ 9t{x,y)dFb^4,{atx,cty) = / 9t{x,y)dH{x,y) + 0{A{t)) 
Je Je 

and in case of (ii) 

/ 9t{x,y)dFb„d,{atx,cty) = / gt{x,y)dH{x,y) 
Je Je 

+ A{t) [ gt{x,y)d^{x,y)+o{Ait)). 



dxdy 



Proof. As f 00, we have to prove that 

^{atx,cty) - h{a 

gtyx,y)^ 

'(o,< 



and by dominated convergence theorem, it is sufficient to show that 
fbt,dt{0'tx,cty) - h{x,y) 



9t{x,y) 



A{t) 



-i>{x,y) 



<K{x,y). 



where K{x, y) is an integrable function. 



□ 
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Let gt{x,y) be the score functions from the pairwise composite hkehhood 
of a max-stable process. Note that J gt{x,y)dH{x,y) = since gt{x,y) is the 
score function. Limit distribution of estimator for dependence parameter can be 
determined by Condition, and Proposition [T]imphes that condition (jl9p should 
be satisfied for the limit behavior. 

We end this section with an example to demonstrate how the proposition 
works. Here we focus on the example with a certain type of gt{x,y), the score 
function obtained from the composite hkehhood of Brown-Resnick process, and 
we intend to show that the condition ([T^ holds assuming that is a bivariate 
normal distribution. 

Example 1. (bivariate normal distribution) Suppose that {X, Y) are i.i.d. from 
a bivariate normal distribution F with mean 0, variance 1 and correlation coef- 
ficient p. First we can prove that bivariate normal distribution satisfies by 
A{t) = (A{t) as t oo). 

Define A{t) — 2ilgt '^"'^ 4'{x,y) — — 2(i-pi) satisfy the condition lil^) . 

Now suppose that gt{x,y) = ^\og fDA{x,y;9) where foA = ^ ^oxdy'^^ ■ ^'^V 
max-stable process can be fitted for modeling annual maxima of data and we 
obtain the score function by our threshold method with the composite likelihood 
approach. Here we arbitrarily choose the Brown-Resnick process with Gumbel 
margins to obtain the joint bivariate distribution of annual data, Fam, cind a 
joint bivariate distribution of daily data, Fda{x, y), is determined by the relation 
ilS^) . With some calculations, the following boundness of the product in 119\) is 
of interest: 

^ x^ -^y^ - 2pxy \ _ 1 \ 
'* 2(l-p2) ; (l + p)2/' 

(20) 

and we can show that \20(l is bounded by an integrable function. See the details 
of proof in Appendix A. 

3.2 Spatial structure and sampling design 

Asymptotic resuhs have been proved for spatial processes which are observed 
at finitely many locations in the sampling region. Central Limit Theorems for 
spatial data have been studied on infill domain and increasing domain struc- 
ture under two types of sampling designs, a class of fixed (regular) lattice and 
stochastic (irregular) designs, in existing literature. Infill domain structure as- 
sumes that the sampling region is bounded and locations of data fill in increas- 
ingly and densely, while the sampling region is unbounded in the increasing 
domain structure. |15j is concerned with more complex spatial structure, called 
mixed asymptotic structure, as a mixture of infill- and increasing domain as- 
sumption. In the mixed asymptotic structure, the sampling region is unbounded 



< 



gt{x,y)exp 



■ exp 
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and sites fill in densely over the region. Covariance parameters are not always 
consistently estimable if the spatial domain is bounded ([30]), while the same 
parameters are estimable under the increasing domain structure (^[18)). Here we 
focus on the increasing-domain case under stochastic design based on setting 
and conditions in [TS]. Increasing domain structure takes advantage of dealing 
with asymptotic properties of estimators easily rather than the infill asymptotic 
structure. We could take account of more realistic setting under the stochastic 
sampling design than the fixed lattice design. 

Suppose that the stationary random field {Z{s); s E M''} is observed at many 
stations s in the sampling region i?„. Under the increasing domain structure, 
Rn is unbounded with n and there is a minimum distance separating any two 
sites for all n. We assume that the sampling region i?„ is inflated by the factor 
A,i from the set i?o, i.e., 

For the stochastic designs of sampling sites, we assume that the sampling sites 
{si, • • • , Sn} are obtained from a random vectors {xi, • • • , Xn} by 

St = XnXi, I <i <n 

where Xi is a sequence of i.i.d. random vectors from a continuous probabil- 
ity density function f{x) and its realization {xi,--- ,cc„} are in Rq. In this 
stochastic design, the sample size n is determined by the growth rate A„ by the 
relation n ~ CAJ^. 

We now consider our threshold approach. Note that we assume the marginal 
GEV parameter rj is known as the simplest case, though we would like to address 
the case t] unknown as well. Assuming that rj is known as t/q, we can rewrite 
(|13p and partial derivatives with the temporal domain fixed, as 

n-l n T 

^ (^) = E E E log ^ (^*' ' ^) 

i=l t=l 
i<j 



de de ' 

d^m _ V \ d^Uj{e) dLt,{e) ( dL.,,{e) V\ 
dede^ ^ Ll{e)\ deoe^ ^ ae \ de 1 j 

where Wij is the weight function on the (i, j)th pair which does not take any val- 
ues outside i?„ , Lij — Fij{xi,Xj)I^Xi>u,Xj>u}~^Pi{xi,u)I^x.^y^,j.^<_y^-^+Fj{u,Xj)I^Xi<u,Xj>u}~^ 

-Fb/iCw, u)/{2;;<t(.2;j<t(} , and Fij = qJ^q^. ■ Here u is the threshold, not a fixed 
constant, which varies as the sample size goes to infinity. 
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We concentrate on the first term of which is the case that both exceed 
the threshold. Let us define notations related with the first term by 

K 
i<j 

dQK{0) _ ^ mjWr 

do f^^F,j{e)' do ^{-•>".-^>"} 
dOdeT f^^ F,,{e)\ dOdQT F,\ 00 J\ de J j^{-.>".-.>"} 

where K is the number of all combination of pairs. 

Next we denote the form of the strong mixing assumption to deal with 
dependence through pairs. Let X{si) = Z{si)I{Z{si) > Un) and J^{Gk) be a- 
field generated by { (X{si),X{sj)); Si, G t/fc, 1 < «, j < n, fc = 1, • • • , K]. For 
any two subsets A and B of M'', the mixing condition is defined by 

H5u52) - sup{|P(A n B) - P{A)P{B)\ : A e F{gi),B e ^-(^2)} 

and let 

d{gi,g2)=ini{\s- s'\ -.seGus'e ^2} 

which is the minimum distance from element of a pair Qi to element of another 
pair Q2- Then the strong mixing coefficient is defined as 

a(a, b) = sup{a(c;i, ^^2) : ^2) > a, Gi,Q2^ T^ai^)} 

where 7?.3(6) = {uf^^-^i : J2'i=i l^il — the the collection of aU disjoint unions 
of three cubes Di, D2 and D3 in Mf^, and it specifies the general form of the sets 
Gi and G2 that are bounded. Assume that there exist a nonincreasing function 
ai(-) such that Imia^oo cti{a) = and a nondecreasing fimction satisfying 

a{a, b) < ai(a)/3(6), a,b> 0. 

In our approach, what we are interested in is the bivariatc function, 

d 

— log F,j{xi,Xj] 9)1 {x, > u,Xj > u) = gk{X{s,),X{sj)) = (21) 

where Zk is obviously different from the original process Z . Let cr(-) denote the 
auto covariance function of the process Zk such that for all Si, Sj,hi,h2 G M'', 

aih) = Cov{Zk,Zi) =. Cov[gk{X{s,),Xis,)),gi{X{s, + hi),X{s, + /la))]. 

Lets^^ = // wfj(Xn^ij)f{xi,Xj)dxidXj = i;?i;|-(A„Xi) where X,y = {xi,Xj), 
Mk = {sup h e M^''} and 7?,. = to simplify the notation, . 

We will use the following conditions which are similar with (S.1)-(S.5) in [15] 
to prove the asymptotic distribution of process. 
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(A'l) JJ\a{x)\dx<^ 

(A'2) Let i?o be a Borel set satisfying Rq C Ro C Rq and i?g be an open 
connected subset of (—1/2, 1/2]"^. The pdf f{x) is continuous, everywhere 
positive with support ^o, the closure of the set Rq C R'^. 

(A'3) Suppose that Xi are i.i.d. from / over Rq and f(xi,Xj) — f{xi)f{xj). 
The joint pdf f{xi,Xj) e [mf,Mf] where nif and Mf are constants in 
(0,oo). 

/ / ''^fj V'^" (^-i '^J ) ) /(^i ,Xj)dXidXj 

q,he M^d 

/A/t;^ //'i'i3(-^n(a;i,3;3))iDi<j(A„(a;i,Xj) + (0,h))/2(3;i,a;j)d3;ida;3 /.\r ,, • • / 

(A 5) ^ — 7 >■ Q2{h) for alli = p,j ^ 

I I ^ij \^n{Xi^Xj]) f {Xi,Xj)dXidXj 

q,he M'^. 

(A'6) = = OiK") for some a € [0, 1/8) 

(A'7) There exist sequences {Ai„}, {X2n} with {Ai„} > {A2„} > log{A„} such 
that 

(i) 7i\(logn)2[^ + A2^] ^o(l) 

(ii) ^f^ilogn^i^) k''-'a,ik) = 0(1) 

(iii) ^ai(A2„)/3(AfJ-o(l) 

(iv) li[XUi{X2n) + EZm„ k'^-'a,{kmXt) o{l) 

Theorem 1. Assume that conditions (A' 1)-(A' 7) hold. Suppose that Zk{s^) 
in \21\) is a stationary stochastic process such that E\Zk{0)\'^'^^ < oo and 
J ai{t)"^ dt < oo for some S > 0. If n/\f^ Ci e (0, oo) as n ^ oo, 
then 

K 

{KslK)-'I^Y.'"K(sk)Zk{s'^) 



fe=i 

d 



N\ 



,a{0) + CiJ a{{0,h))Q2{h)dh + Cf J a{h)Qi{h)dh\ 



Proof. Proof of Theorem [T] is shown in Appendix C. □ 

4 Asymptotic properties of spatial dependence 
parameter estimates 

4.1 Asymptotic normality and consistency 

We use the following regularity conditions to obtain an asymptotic behavior of 
estimates of dependence parameters. 
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(Al) The support x of the bivariate density function of the data does not 
depend on e 8 and the parameter space 8 is an open subset of M.P with 
identifiable parametrization. 

(A2) The pairwise composite log likelihood is at least twice continuously differ- 
entiable in 6. 

(A3) (smoothness of composite likelihood) Qt{S) exists and is continuous and 
M{9) is also continuous in a neighborhood 8* of Oq. 

(A4) For all Oq E Q, there exists an integrable function M{x, y) such that 
d^QT{e-,x,y) 



sup 



86,86. 



< M{x,y), ij = I,-- - ,p. 



(A5) The third partial derivatives of the composite likelihood are bounded by 
integrable functions. 

(A6) (equivalent condition of Proposition [T]) The score function of composite 
likelihood D satisfies that 



Dix,y) 



fbt,dAatx,cty) - h{x,y) 



Ait) 



ipix,y) 



<K{x,y) 



which K{x, y) is integrable. 



Theorem 2. (Asymptotic Normality) Suppose that condition Hl^) with (i) or 
(ii) is satisfied and conditions of Theorem]^ hold. Suppose N — >■ oo, {bk,dk) = 



{h{k)N,d{k)N) {xo,yo), and A{kN) = 0[ 



^JNsl^A{kN) \e [0,oo), 



and either A = and (i) holds, then the the solutions of likelihood equations 
verify 

Vn{sIn)-'^H§ - Oo) A A^(O,H(0o)"M^o)H(0o)"'), (22) 
or (ii) holds, then 



NislN)-'/^ie-eo) A 7V(H(0o)^'b,H(0o)^'V(0o)H(0o)"'), (23) 

where H(0o) = IE[-D'(0o)], b = lim^^oo (^s^j^j^'^'Ej Ef=i D(^o; ^^'^)} 
(defined below) and V(6>o) = E[D(6»o)D(0o)'^]- 
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Proof. Denote that 

K 

K 

+ Fj{u, 2^j)-^{a;i<u,a;3>u} 
K 

:=^^Wij\og{Liij{d;Xi,Xj) + L2ij{d;xi) + L3ij{d;Xj) + L4ij{d)}, 



dO ~ LiAO)' dO 

K 

d 

— {Luj{e;xi,Xj) + L2rj{6]Xi) + L3^j{d;xj) + L4ii(^)} 

N 

.(fe) yik). 



k=l 
N 

J2 {Di{e) + D2{0) + D3{d) + D4e)} 



k=l 



where 



i<ij t J \ / 

and K is the number of all combination of pairs. Wc now consider the 
number of exceedances as a primary role in deriving the asymptotic behavior. 



17 



By notations and condition (A2), we have Taylor expansion about 6q as follows. 
dQrie) 



- 



80 

dQrie) 



de 



+ {e- 00 



902 



(0-00) 



where 9* lies between and Oo 



J2 D(0o; \ )) + (0 - Oof D'(eo; 

k = l fc = l 

+ ke-0ofJ2^"{0*; ))(0 - 0o). 



Then we rewrite the equation as 

N 



fc=i 

I k=l k=l J 



and then 



iED'(^o)4(e"-^„fif:D"(r) 



/c=l 



fe=l 



V^^SlAT fc=l 



(a) 



(b) 



(24) 



We establish the following for separate terms in equation (|24)) : 



(I) By the consistency of 9 and condition (A5), expectation of the last term 
in parentheses can be ignored. Since 9 is consistent, 9 E Q* with Peg- 
probability 1. Let Be 8* be a closed ball with the center 9o. By the 
condition (A4), 



sup 

SeB 



d^QT{9o-x,v) d^QT{9-x,y) 



is bounded and then, for large N , 

N 

N 



lim sup^ 



N 



^Ed'(^o)-^Ed'w 



k=l 



fc=l 



< £ 
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in probability (see details in [12] )■ jf J2k=i ^'(^o) — D'(0o)] converges 
to by the law of large numbers, and hence (a) converges to E[— D'(0o, Vo)] 
in probability. 

(II) First consider that (X^^\ X^j'^), k — I, ■ ■ ■ , N are i.i.d. from exact mul- 
tivariate GPD distribution H. 

N 



N 

(no bias). 



Then by Theorem[l] (b) converges in distribution to N{0,Y{6a)) where 

V(0o) = E[D(0o)D(0o)^] 

= E[(Di + D2 + D3 + D4)(Di + D2 + D3 + D4)^] 



Var[DiD/] + VarpaDa^] + VarpsDg^] + Var[D4D, 
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Var[DiDi^] =a(0o;O) 

+ Ci J a{eo: {0,h))Q2{h)dh + Cf J a{eo;h)Qi{h)dh, 

and VarpaDa^], Var[D3D3"^] and Var[D4D4^] have similar forms with 
the variance of DiDi"^. Note that the event {xi > u,Xj > u} of Di is 
uncorrelated with the event {xi > u,Xj < u} of D2, and Cov{Di, Dj) = 
for i ^ j. 

(k) (k)\ 

Now suppose that [X^^ , X^j ), k — ,N are from Fbi^^dk not H. 

U F £ D{G), there exists the exceedance level {bk,dk) such that Fb^,dk 
converges to H as {bk,dk) — >■ {xQ,yo)- The bivariate generahzed pareto 
distribution H preserves under the suitable change of exceedance levels 
(see [53]). 

The second-order condition (fTSl) describes the difference between F^^^dk 
and H with the remainder function A{k), i.e., as fc — )■ 00, with the second 
order condition (ii) 

limfiup\Fb^^dk{'^kXi,CkXj) - H{xi,Xj)\ = 0{A{k)). 

Proposition [1] (ii) results from the condition (A6), and by the property of 
score function 



1= /ED(0o;XW)dF,,,,,( 



0V^A(fc)--V / 5]D(0o;XW)d^'(x„x,) + o(A(fc))^A/z, 
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where fi = liniAr^oo / D(0o; XW)d^'(a:,, x^). 
Then for some finite vector b 



fc=l 



and (b) converges in distribution to N{b,Y{6Q)). Therefore the limit dis- 
tribution of 0, follows by Slutsky's Theorem. If the second-order 
condition (i) holds and ^/Wsf^A{kN) — > 0, 6 = which implies no bias 
and then (Ull) holds. 



□ 

To prove consistency, we describe the theorem of [1]. 
Theorem 3. (Amemiya, 1985) Assume the following: 

(Bl) Q is an open subset of Euclidean p- space (the true value 9q is an interior 
point of Q ), 

(B2) The criterion function Sn{0) is a measurable function for all Oq G Q, and 
VS'at exists and is continuous in an open neighborbood of 9^, 

(B3) j^Sn{6) converges in probability uniformly to a non- stochastic function 
S{9) in an open neighborhood of 9q, and S{9) attains a strict local maxi- 
mum at 9o ■ 

Then there exists a sequence cat ^> such that 

P{39* such that \9* - 9o\ < ew, VS'Ar(r) = 0} ^ 1, as N ^ oo. 

Theorem 4. (Consistency) Let xl*"^ = (X^'f \ Xj^j"^), k ^ 1, - ■ ■ ,N be i.i.d. 

random variables with bivariate distribution F. Let 6 be the maximum pairwise 
composite log-likelihood estimator such that 

K T „ 
o 



i<'J t=l 



9=e 



If the second moment condition of composite score function is satisfied and con- 
ditions (Al), (Bl) and (B2) hold, then there exists 6 such that \0 — Oq\ < tN 
and VSNiO) — for any sequence cn ^ 0, as N ^ oo and {bk,dk) = 
{b{k)N,d{k)N) {xo,yo)- 

Proof. Assumptions (Bl) and (B2) in Theorem [3] are satisfied by our criterion 
functions and assumptions. Jensen's inequality implies 



20 



We rewrite it as 



Ea 



log 



< -i^ 9q ^ arg max Eg,, 
see 



lOE 



fix; 9) 



f{x;9o) 



Here a sum of pairwise log-likelihoods can be considered. Let 

K T 

Sn{0) = 5^5]w.,logL(Xt„Xt,;0). 

Kj t=l 

We know that by the law of large numbers, 

N N 

j^SNiOo) = ^5:z.If logL(xf,x(f ;0o) = j^Y.^^'nogL{X^,'^ ■,60) 

k=l k=l 

-^Ee„{w^^hogL{X^,^^;eo)) =: 5(0o). 

By the moment condition of V5Ar(0), we have that E\V Sn{9*)\'^ < Cq for some 
Cq. Using a Taylor's expansion, 

(l^Ar(0)-5(0)) - {^SNido) - SiOo)) 



< ( -E\VSN{9*r 



E\vsN{e* 



\9-en\\^ 



-^CoWe-eoW" (26) 

for some 6* and 9** between 60 and 6. By the moment condition of VS'Ar(0), the 
right hand side of ([26l) converges to uniformly over a sequence of 1 1 — 0o 1 1 < 
cn as Eat -> 0. Also we have that j^SN{6a) — S{9q) by the law of 
large numbers and -^S'Ar(0) converges in probability uniformly to S{6) on a 
neighborhood of ^o- 

Now we claim that S{6) attains a local maximum a.t 6 = Oq. The previous 
result (PSI) implies that 



E, 



log 



Y{kL{x[^^;e) 



YlkLixt^-.e, 



< iog Eg, 



UkLixl'^;0) 
UkHx^t'^^Oo) 



= 



and for any 6, 

£;«„(logni(^f ^0o)) > Eeo{\ogllL{X^,'^;e)). 

k k 

where the equality holds with (Al), the identifiability assumption of parameter. 
Eg,, \ log l^*" '^'•^A.^ 1 < holds for any distribution of Jci'^-' with finite second 

moments of score function, and the maximum of -E'eo[log]^j. L{x\^'' ] 6)] over 6 
is attained at = 0o- Thus we prove the (B3) of the Theorem |31 □ 
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4.2 Simulation 



We conduct some simulation studies to illustrate the asymptotic behavior of the 
estimators described in Section 14.11 The simulation is examined for the daily 
max-stable process with unit Frechet margins with T — 1000 days during 10 
years, i.e., M — 100 in equation p^ . We consider the Gaussian extreme value 
processes with two different spatial dependence structures of the covariance 
matrix: 



(i) the Gaussian extreme value process with Ei (a = 2, /3 = and 7 = 3); 

(ii) the Gaussian extreme value process with E2 (a = 2, /3 = 1.5 and 7 = 3). 

We generate n — 20 stations from the uniform density function /(•) over Rq = 
(—1/2, 1/2]'' and determine the growth rate A„ — s/n in case of c? = 2 to satisfy 
the relation n ^ in the spatial structure and stochastic sampling design 
of sites. To adjust the threshold approach based on the pair wise composite 
likelihood, we consider a weight function such that for some constant i5o, 



where /i is a distance between two stations. Here (Sq is selected by \/2n/2, the 
half diagonal of sampling region, which satisfies the condition (A'6) on growth 
rate of weight function for the asymptotic result. 

To illustrate the asymptotic performance of estimates for dependence param- 
eter Q = (a, /3, 7), the averages of the estimators are compared to the asymptotic 
mean of 0. In each model, the estimation of dependence parameters is based on 
500 replications, and the classical Monte Carlo integration is used to implement 
the theoretical bias and variance of the estimators as the number of exceedances 
N increases. 

Theoretical bias and average bias of estimators d for Smith model (i) are 
plotted in Figured] As the number of exceedances increases, bias of estimators 
(gray curve) tends to decrease towards the theoretical bias (solid curve) though 
each estimator shows the different slope on the decay. The bias of a goes on 
with the pattern of decay of theoretical one, while bias of ^ and 7 decreases as 
theoretical bias goes up to the line of zero bias. 

This irregular pattern of each dependence parameter estimation might be 
caused by the interaction between parameters in estimating them as components 
of covariance matrix. Now we plot the extremal coefficient curves with the 
parameter estimators and compare them with those estimated directly. One can 
expect the problem to be reduced when working with the extremal coefficient. 

Figure [2] shows estimated extremal coefficient functions by d. As the num- 
ber of exceedances increases, the color changes from cyan to magenta. Extremal 
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Figure 1: Graphical summary of asymptotic behaviors of d, /3 and 7 for Smith 
model (i) from left to right. Gray curve is the average bias of estimators, gray 
dashed curves are the boundary of 95% confidence interval, and black solid curve 
is the theoretical bias. 




Colors Cyan(')-Ma9enta(1U00) 



Figure 2: Extremal coefficient functions for the Smith model (i). Upper thin 
color layer is based on theoretical mean of estimates and lower thick color layer is 
based on average estimates. In a layer, each line represent a extremal coefficient 
curve at each N and the line changes the color from cyan (TV = 1) to magenta 
{N = 1000). Black solid line is the true extremal coefficient curve. 
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Figure 3: Graphical summary of asymptotic behaviors of d, /3 and 7 for Smith 
model (ii) from left to right. Gray curve is the average bias of estimators, gray 
dashed curves are the boundary of 95% confidence interval, and black solid curve 
is the theoretical bias. 




Colors Cyan(')-Ma9enta(1U00) 



Figure 4: Extremal coefficient functions for the Smith model (ii). Upper thin 
color layer is based on theoretical mean of estimates and lower thick color layer is 
based on average estimates. In a layer, each line represent a extremal coefficient 
curve at each N and the line changes the color from cyan (TV = 1) to magenta 
{N = 1000). Black solid line is the true extremal coefficient curve. 
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Figure 5: Mean squared error of d, /? and 7 for Smith (i) from left to right 



coefficient by the asymptotic bias overlapped almost with the true coefficient 
function (black solid curve). As the number of exceedances increases, The ex- 
tremal coefficient curve measured by dependence estimators approximates the 
theoretical extremal coefficient curve. However, there still exists a gap between 
the theoretical extremal coefficient and estimated one and the gap gets broader 
as the distance between two locations is larger. 

Theoretical bias and average bias of estimators 9 for Smith model (ii) are 
shown in Figure [3] As the number of exceedances increases, bias of estimates 
tends to go towards the pattern of theoretical bias. There is some gaps be- 
tween theoretical bias and estimated bias though the estimation of dependence 
parameter is much more stable comparing with that in model (i). 

Figure |4] shows estimated extremal coefficient functions by 0. As the num- 
ber of exceedances increases, The extremal coefficient curve measured by depen- 
dence estimators approximates the theoretical extremal coefficient curve. Unlike 
the gap in Figure [3l the estimated extremal coefficient is catching up with the 
theoretical one along by a little gap. However, the quality of asymptotic approx- 
imation seems dependent on the degree of correlation /3 since Figure [5] shows 
the poor approximation to the true extremal coefficient curve comparing with 
Figure [H 

Suggestion on the choice of the threshold point is discussed further now. For 
the simplicity, the threshold can be selected as the value of the 95th percentile 
of distribution function in practice. However finding an optimal threshold is 
another important issue and we suggest an optimal threshold minimizing the 
mean squared error, which incorporates both the bias of the estimator and its 
variance based on the asymptotic normality in Section [4.1l 

Figure [5] shows the mean squared error for each estimator in Smith model (i). 
The mean squared errors of d, /?, and 7 are decreasing rapidly against TV and 
show the stability between TV ~ 3500 and N = 5000. As shown in Figure [2l the 
theoretical extremal coefficient has a nice approximation to the true coefficient 
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Figure 6: Mean squared error of a, (3 and 7 for Smith (ii) from left to right 



function, and the increases of squared bias seem to be less effective than variance 
decreases on the selection of threshold to minimize the MSE. 

The mean squared error for each estimator of Smith model (ii) is shown in 
Figure [6l The mean squared errors of a, /3, and 7 are decreasing rapidly as N 
increases to 1000, and have the minimum between N = 1500 and N = 2000. 
In Figure SI the theoretical extremal coefficient shows the poor approximation 
to the true coefficient function. Thus calculation of MSE is affected by the 
increase of bias as number of exceedances becomes greater than 1500. The 
threshold point is suggested as the value between 9Qth and 95th percentile. 

5 Discussion 

The threshold approach takes advantage of avoiding the loss of information 
which is caused when we are concerned with only maxima of data. Our method is 
expected to become one promising tool to characterize the dependence structure 
in spatial extremes, we have suggested the modeling of the bivariate exceedances 
over threshold and it leads to a simplified dependence structure for max-stable 
processes. An important motivation of this methodology is the possibility of 
threshold approach to construct approximation of the joint distribution, by 
assuming an asymptotic distribution of exceedances over a given threshold. We 
have derived our simulation results under two Smith models to examine the 
asymptotic property of estimates. 

Moreover, we have also investigated an optimal threshold to minimize the 
mean squared error based on the asymptotic behavior of the estimator for de- 
pendence parameter. The choice of optimal threshold would be an open topic 
itself for further research. It provides very valuable information in the field 
of environmental statistics. When we are interested in flooding, for example, 
which may be considered as extreme events, choosing the adequate threshold to 
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avoid the risk of flooding might be useful for quantifying the spatial extremal 
dependence. 



Appendix 

A. Example of Proposition 1 

Suppose that {X, Y) are i.i.d. from a bivariate normal distribution F with 

mean 0, variance 1 and correlation coefficient p. First we would like to prove 
that bivariate normal distribution satisfies (18) in the paper. We consider G 
in (14) as a bivariate extreme value distribution with Gumbel margins, and 
suppose the limiting form of bivariate normal G{x,y) ~ exp{— — e~^} in 
the case of the independence. A max-stable process with unit Frechet margins 
will be fitted and the transformations X' = log X and Y' = log Y can be made 
from unit Frechet to Gumbel. 

Mills ratio for a normal density implies that 

1 - $(.t) r 1 1 1-3 1 • 3 • 5 
(j){x) \x x^ x^ x"^ 
P{X>x,Y>y) (l-p2)' ^ 



(l){x,y) {x- py){y - px) 



{ 



I-(I-P^) 



1 P 1 



(.T - py)2 (x - py){y - px) {y - px)' 
(see [?] for the bivariate normal density). From the fact that 

1 - F(x, y) = \- ^(x) + 1 - $(y) - P{X > x,Y > y), 
we could set the lower bound and upper bound for ^"^"^'^^^ such that 

l-F{x,y)Y ^ l-F{x,y) ^ (l-F{x,y)^'' 



where 



<l^{x,y) ) ^{x-.y) V <\>{T'^y) 

l-F(x,2/)\^ 11 1 1 {^-P"? 



(j){x,y) J X y x^ y^ {x ~ py){y - px) 
l-F{x,y)Y ^ 1 (1-p^)^ ^ (l-p2)3 



<P{x,y) J X y {x-py){y-px) {x - py){y - px) 

^ P + ' 

(x-py)'^ (x-py){y-px) {y - px)'^ 

From the well-known results of extreme value theory, define 6t by 1 — $(6t) = 
i and at = ^/bt- Or we might set normalized constants 

1 

at 



VTlogi 

I. An — I I (log log * + log 

bt = V2^t ^_ . 
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Conditional distribution of exceedances over threshold is written as 

^ , . , t{l- F{atx + bt,cty + dt)} 

Fb,,dMtx,cty) = 1 -J- . 

t{l - F{bt,dt)} 

and we now concentrate on i- J'(«tx+bt,cti'+'^*) ^ 

1 - F{atx + bt,cty + dt) 1 - F{atx + bt,Cty + dt) 4>{bt,dt) (t){atx + bt,cty + dt) 



> 



1-F{bt,dt) (l){atx + bt,cty + dt) 1-F{bt,dt) (l>{bt,dt) 



x + by y + dl {x + bjf {y + blY + + 
&t \ (l){atx + bt,cty + dt) 



j l263-(l + p)26|+(i±^iM/ <j>{btM) 

and also, 

l-i^(ata; + 6(,C(y + dt) 



< 



l-F{budt) 

' h dt (1-P')'6| 



l-(l-p2)x 



{x- fyy + bl{l- p)) {y- px + b}{l- p)) 
'I 



b't 



2bl - (1 + pfbt - 2 



{ 



{x-py + 62(1 - p)) {y-px + bl{l - p)) 
^{atx + bt,cty + dt) 
4>{h-,dt) 

26? - (1 + + M^to) W ^,3 . ^(^/J,^ + f,^^ yl^^ + 



Thus 



b\ ' /l26|-(l + p)26t-2j ^{btM) 



T? I \ uf \ 1- F{atx + bt,cty + dt) , . _x , -y. 
Fb^,dAatx,cty)-H{x,y)= 1 - Fibt dt) + e 

f 26? - (1 + p)26^ - 26, U <^(x/6t + 6,, y/6« + bt) \ 

~l 263-(l + p)^^? + ^i±€f^ A <i>{bt,bt) ^ ^ i 

and 



26?-(l + p)262-26t 



26 



(l+p)"(2-p) 



26? - (1 + p)26| + ii±^ie^ 26? - (1 + p)26f + (l±^i^ 
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We obtain the formation of (18), 

lim P^.MatX,Cty) -H{x,y) ^ 

t-foo A(t) 

where Ait) = ^ = and ^{x, y) = exp { - f±^} + e"- 

Next, 

„ . ^ atct 1 

Jbt,dt{atx,cty) 



exp 



1-F{bt,dt) I'K.JT^ 

(atx + htf + (cty + dtf - 2p{atx + bt)icty + dt) \ 
2(1 / 



(t){bt,dt) (j){atx + bt,cty + dt) . 4>{bt,dt) , 

= F/L ^ A . /, , N = atCt- ■ Vt{x,y) 

1-F{bt,dt) (t>(bt,dt) 1-F{bt,dt) 

(27) 

where (j^ix, y) is a bivariate normal density with correlation p. 

a factor of fbt,dt (itx, Cty) in the equation \21\ has the lower and 
upper bounds that 

4){x,y) \^ ^ Hx,y) ^ f 4>{x,y) 



l-F{x,y)J - l-F{x,y) - \l~F{x,y), 
where 

(p{x,y) V _{l , 1 (l-p2)2 (l_p2)3 



l-i^(a;,y)/ [a; y {x - py){y - px) {x - py){y ~ px) 

1 p 1 ^^-^ 



[x-pyY {x-py){y-px) {y - pxf 

Hx,y) f 1 ^ 1 1 1 {1-pY 

l-F{x,y)J \x y x^ y^ [x ~ py){y - px) 

Since fb^^dt{o.tX, cty) = atCtJz^FJ0d^^^t{x^ y), using above normalized constants 
and assuming bt — dt 



atCt 



atCt 



Hbt,dt) \^ bf 



l-F{budt)J 2&f ~(l + p)262+ (i±ig^ 



(j){bt,dt) \^ bt 



l-F{bt,dt)J 2&2_(l + p)26,_2 

(j>{x/bt + bt,y/bt + bt) ( x'^ + y'^-2pxy x + y' 

Vt{x, y) = = exp ■■' 



mM) "I 2(l-p2)52 i + 
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Thus we could get the following form of bounds 

rL , N_ ^1 0(x/6t + &/,.?y/??t + &t) 

h,A. ^<^tx, cty) _ ^ ^^^^^^ ^ ii+p)H2-p) ^(6„ 6,) 

fU ( ■s_ 4>{x/bt + bt,y/bt + bt) 

f»uiMtX,cty) - _ (1 +p)2^^ _ 2 ^(btM) 

Meanwhile 

Therefore 

fh,dt{atx,cty) - h{x,y) > {fbt,dt{atx,cty) - h{x,y)}^ 

^ &t 4>{x/bt + bt,y/bt + bt) 

263 _ (1 + ^)252 + (i+pP(2-p) </,(6„ 6,) 

fbt,dt{atx,cty) - h{x,y) < {fbt,dt{atx,cty) - h{x,y)}^ 

^ bt 4>{x/bt + bt,y/bt + bt) 

262-(l + p)26t-2 cl>{bt,bt) 

Define A{t) = 2ilgt (^(*) ^ as f — )• 00) and il){x,y) = — 2(1-^2) to satisfy the 
condition (18). Then we could show that 

fbt,dt{atx,cty) - h{x,y) ^ ^ fbt,dt {atx, cty)^ - h{x, y) , ^ 

A{t) n>yx,y)^ 1/(2 logt) ^^""'^^ 

/ x + 2/\r6t / ^x^ + y'^ - 2pxy\ 

~ '"H " TT7 J 1 2 r 2(1 -p^) j 



2(1 J (1+^)2 

fbudt{atx,cty)-h{x,y) ^ ^ hudMtx, cty)^ - h{x, y) ^ ^ 

A{t) nx,y)S 1/(2 logt) '^^'"'y' 

( x + y\(b, ( 2X^+y^-2pxy\ 1 

This limit for bounds of '^'''^'(^^ ^' — "0(2^: J/) will be used to prove that the product 
of a function gt{x,y) and '^'''^'^'^-j ^ ~ ^(2;,?/) is bounded by an intcgrable func- 
tion as shown in (19), Proposition 1. Suppose that gt{x, u) = -§0 log foAix, y; 9) 

where foA = ^ ^dxdy'^^ - max-stablc process can be fitted for modeling 

annual maxima of data and we can obtain the score function by our thresh- 
old method with the composite likelihood approach. We arbitrarily choose the 
Brown-Resnick process with Gumbel margins to obtain the joint bivariate dis- 
tribution of annual data, Fam-, and a joint bivariate distribution of daily data, 
FDA{x,y), is determined by the relation (12). 

Fam{x, y; 9) = exp{B(a;, y; 0)}, 



30 



wherei3(.,y;^?) = {-l$(^^p^^ 



and 



where J{x, y; 9) 



{h-0) ^^J y \ 2 ^^KM) yJ) 



log fDA{x,y]0) = —B{x,y;e) + log J{x,y; 



_ 1 d^B(x,y-e) 1 dB{x,y;0) dB{x,y-e) 



Therefore, 



9t{x,y) 



M dxdy ~^ dx ' dy 

1 dB{x,y-e) , ^,,^,„i/9J(a;,y;0) 



where = 



M ae \ 



M de 



d^B{x,y;e) 
dxdy 



+ j{x,y-ey 



80 



1 d_fi 
TT^ de \ 



dB{x,y;e)\ dB{x,y\6 



dx 



dy 



(28) 



1 dBjx.yfi) 
JI^ dx 



d_ ^B(x,y,e) ^ ^ With some calculations, the derivatives of J{x, y; 9) and B{x, y; 9), 

shortly J and B, can be obtained as in Appendix B and the boundness of the 
product is of interest: 



9tix,y) 



fbt,dt{atx,cty) - h{x,y) 
A{t) 

< 9tix,y)exp ( - 



1+pJ [2 



tp{x,y) 
bt 

■ exp 



1 



^ x^ + y^ - 2pxy \ 

2(1 ; (i+p)2 

(29) 



Case (i) x ~ y: 

dB 
'89 



89 



1 



M 



Ji9) = y^e-^^{^]+j^e 



2^7 
1 



V2 

2 



'09 



89 



1/1 



V2 

2 

1 / 1 



Then 

9t{x,y) < 



8-f 
'89 



M2 ^2^ 



-2x, 



2 



and therefore, for some constants Ci 
ht,dt{atx,cty) ~ h{x,y) 



9tix,y) 



<Ci 

<C2 



ipix,y) 



e i+p < bt exp 



A{t) 

-X 

— - 1 p ^ p 1+p 
89 , 

\/2x 



(\ + p)hi) (i + p)2 



btii + p) 



.V^Tpbt ^2(1 + p), 
which imphes that ([29| is bounded by an integrable function. 
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Case (ii) y = x + k and a: — > oo: 
Let 



2 

1 1 



4^7 2V/73 



4^7 

dB__ /dj 
'80 " \de 



{x - y) 
{x - y) 
{y-x) 



V7 k 

2 V7 

27 V 2 V7 

J_/%/7 

27 V 2 V7 
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^(^i-^(60(a)+e-'=a0(6))|, 



1 



a + e 



0(a) 



0(a) 
V7 



V7 
V7 



$(&) 



V7 
0(6) 
V7 



'de ^ I 9^ 



— e--(0(a)fc3(fc) + e-V(&)fc3(-fc)) 



+ ^{0(«)fci(fc) + e-^mk2{-k)]^^e-^l^m + ^) 



0(6) \ 0(a) 



-2a; 



M2 



■{0(a)fc2(fc) + e-"0(6)fci(-fc)} $(a) + 



0(fl) 



V7 



V7 / V7 

V7 J. 



where fci, fc2 and fca are defined in Appendix A. Then for some constants 

9t{x. y) < (^g^ I _ _L (^^ j (Xi60(a) + i^2a0(6)e-'^) je^ 

and therefore, for some constants Ci 
' htMt{atX,Cty) - h{x,y) 



9t{x,y) 



A{t) 



< C*! ( I e e i+p • - exp <i - 



tp{x,y) 

2(1 - p)x'^ + 2(1 - p)kx 



< ^20 



89^ 

f 2x+{3 + p)bj + k 



2(l-p2)52 



V ^2{\ + p)bt 
which implies that (1291) is bounded by an integrable function. 
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For the general case of a; ~> oo and ?/ — > oo, the boundness can be obtained 



In the first term consists of the components; — e ''^(py^ + ^^-^ 

and -e-y^[^ + ^ j. In the second term of gt{x,y), J(a;, y; 6*)-! 



is also dominated by (T'^^y^ + and e"?'?!)^^ + ^j. Then ([291) is 

bounded by a function of (/>(Cix, C2y) for a constant Ci, which is integrable. 

B. Derivatives of the Functions 

dB _ ( dj\ f p-x^f , y - ^\ f 1 y ^ ^ 



V7 / V4^/7 2v/t 



9a; V2 /7y^V2 ^/7/^/7 



d0\dx \de 





a; - y 








y- a; 










27, 





^ V7 / V8V7 2V^ 27^^ 



e 



where k\{x) 
and fc2(a;) 



dd)\ "V 2 77/ ' ' "V 2 77 

1 1 : 



8^7 2Vr 2Vr 27Vr 

1 1 a;2 



8^7 2^7^ 27v9 



a _ ( dj\ \ ^^y.( y/^ ^ X - y\ f 1 1 x-y ^ {x-yf 



de[dy) U^Jr ""^V 2 ' V7;V8V7 2^ 2^' 27v^ 



2 V7 / V8V7 2Vy 27VT 
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dxdy 

d / d^B \ 
de\dxdy) 



97 

'de 



2 V7yV2V7 ^ 



V7 



r 1 i_ f 1 3 



2 ^/7A2^/7 ^ J' 

{y - xf _ {y - x)^ \ 
47-^/7^ 272-^73/ 



iy-x) + 



{- 

97 



1 1 

16Vf ~ 4^ 



+ 



7^ 27 y 7' 



3 ^(^_^) + (^-^/)' i--yr 



47^7^ 27 V7^ 



where ^3(0;) = — 



1 



16^7 4V7 



1 



7^ 27 y 7' 



x + 



47 Vt^ 272-^/7^ 



Leta=^ + ^and6=^ + ^. 

- ^axdy ^ dx' dy 
1 



^,e--<^(a)6 + e-M6)a 



x/7 



%/7 



%/7 



5J _ 1 d/d^B\ 1 d /dB\ dB 1 d /dB\ 

'de ~ Mdeyoxdy) M^Mva^J ' ay M2 ax "de\dy) 



1 



e ^(/>(a)A;3(y - x) + e ^ (t){h)kz{x - y) 



+ -j^^e-^mMy -x) + e-^mUx - ?y)||e-^ (^$(6) + ^ 
+ ^|e-^<^(6)fci(x - y) + e-^4>{a)k^{y - x)| je"- (^$(a) + ^ 



C. Proof of Theorem 1 
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WLOG, assume wk = Wij 

= ^ ^ ^ ^ Wij {Xn{Xi, Xj))Wpq (A„(Xp, Xq))a{\n{Xi, Xj), Xn{Xp, Xg)) 

i j>i p q>p 

= JU2J2Y1 (Xij, Xp,), = {xi,xj) 

i j>i p q>p 

Assume that f{xi,Xj) = f{xi)f{xj) € [m/,M/] where m/ and Mf are con- 
stants. 

/ / Wij{\n{Xi,Xj))Wpq[\n{Xi,Xj) + h)P{Xi,Xj)dXidXj 



< 



J J wfj{Xn{Xi,Xj))f{Xi,Xj)dXidXj 
J J Wij{Xn{Xi,Xj))wpg(Xn{Xi,Xj) + h)dXidXj 
'm-f J J wfj(Xn{Xi,Xj))dXidXj 



f IS w^^{Xn{xi,Xj) + h)dxidxj 



^ — h f f 2 \ : (by C-S inequaUty) 

J \ j j Wfj [Xn{Xi, Xj))dXidXj 



< — i < oo. 



ruf 

Ea%: = K{K - l)EwK{Xn^ij)WK{Xn^pq)cr{Xn{'Xij - Xpg)) 



+ KEwKiXnyi^jf<7{0) 

n(n-l)(n-2)(n-3) 



EwK{Xn'^ij)WK{Xn^pq)o-{Xn{'X.ij - Xp,)) 



4 

+ n(n - l)(n - 2)£;wK(A„Xy)wii:(A„Xig)cr(A„(Xjj - Xj,)) 

+ !^(!LlI)^^^(A„X,,)V(0) 

n(n- l)(n-2)(n-3) 



4 

f{Xij)f{Xij+X-'h)dXijdh 



'^n^'^ J '^(h) j u;y(A„Xy) 



+ n{n - l)(n - 2)A-'' y ct((0, h)) x 
y u;,,(A„Xy)«;,g(A„Xi,- + (0, /i))/(Xi,)/(Xi, + A-i(0, /i))dX,,d/i 



+ !^(!L^^^^(A„X,,)V(0) 



ii'n2A-2''£;M;|f(A„Xi) J a{h)Qi{h)dh 

+ KnX-'^Ewl{Xny^i) j (j{{Q,h))Q2{h)dh + KEwK{Xny^if<j{0) 
--{Ks\^)(^a{Q) + Ci j u{{Q,h))Q2{h)dh + Cl j a{h)Qi{h)dhj , 
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as n -T> oo {K — > oo), by (A'4), (A'5) and dominated convergence theorem. 

K K 



a=l b=l 



hiK{x) = Ehiiix, Xi), X e M^d 
(by Eq. (5.6) in [H]) 



K 



- = ^ [/ik(X„ X,) - £;;i;^(Xi, Xi)] 

a=l 
K-l 



fc=i 

K a-1 

a=2 6=1 
DiK + D2K + D^K 



\EhK{x,XiY\ = 



< {MlfMjX-^'' j \a{s)Yds, by (A'l). 
\EhK{X^,X2Y\<E\EhK(X^,X2r\X^)\ < C{Mf,a{-))M^'-X;,^'' 



Then 



K=l 



Y: E{a\ - Eaj,)'/{K'X-''EwliXM) 

00 

<C{Mf,a{-),C^)J2 

00 

= C{Mj,a{-),C^)Y,h- 



Ml 



'^Vs<(A„Xi)y \K^Xn'"' 

fK^X-^'' 



2 n4 



K=l 



00 / 2 \4 

C(M/,a(.),Ci)^npL 

_R'=l 



< 00, by (A'6) 
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since 



K K 

ED^ < C[Y,{K - hfEhiK{Xif + { Y,{K - afEhiK{X^fY] 

6=1 a=l 

<C{Mi,a{-))KHll\~^\ 

K a-1 



a=2 6=1 

<CK± 



a=2 



E{{a - l)E[{hK{Xa,X^) - /iik(Xi))4|X,] 

+ ((a - l)E[{hK{Xa, Xi) - /ii^f(Xi))2|X,])'} 
+ E[{a- 1)(/iik(X„) - ^/iik(Xi))] ' 

(see details in Eq. (5.7)-(5.9), 15 ). It follows the analogous result by Lemma 
5.2 (i) in [15 . If n/Xi Ci € (0, oo) and (A'l), (A'4) and (A'5) hold, then 



cj{0) + Ci j <j{{0,h))Q2{h)dh + Ci J cj{h)Qi{h)dh^ . 



Let a = = - EZkis''). Define for c > 0, 

Vk^ikI{M<c)~Ei^I[M<c) 

ik = ^ki{M>c)-Eioim>c) 

where = ^fe(O). Let S]^ = Z]f=i Wkrik, Sj^ = Y.k=i Wklk, and 

ai*(a;;c) = Cov{(,kim < c),^oI{\^o\ < c)), 
a;{x;c) = Covi^kim > c),Co/(ieo| > c)). 

We separate the sum of centered processes into two parts, 

K 



SK = J2wkis'')^kis'') 

fc=l 

K K 



fc=i 



fc=i 
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By the moment condition on Zk{-) and the strong mixing condition, 

poo 

Jo 

|CT*(a;;c)|da; < oo Vc>0,j = l,2. (30) 
Since |(5i(a;)| < 1, we obtain that for all x and c, 

J J a{x)Qi{x)dx — J J al{x;c)Qi{x)dx 

< 1 1 l^\Cov{^k{x)I{Mx)\> c),^k{0))\ 

+ \Cov{^k{x)I{\U=^)\ < c),efc(0)/(|a(0)| > c))\jdx 

< C{d){EMO)\'^')^ {E\UO)I{\UO)\ > 

nOO 

Jo 

— > as c — )■ oo 

y" J a{x)Qi{x)dx - J J crl{x;c)Qi{x)dx ^ o{l). (31) 
The similar one can be applied to the form with Q2. 

P ^lini^ [ES],* (c) - ali,{c)] /{KsIj,) = = 1 (32) 

where 

^2,k{c) = K{K - l)EwK{Xn^l)wK{K^2)cr2 (A„(Xi - X2); c) 
+ KEwl{XMa;{0;c). 

From the previous proof of the asymptotic variance, we can obtain the result 
such that for any c > 0, 

^Ik{c) = {KsIj,)I^C! J a*2{h;c)Qi{h)dh + Ci J a;{{0,hy,c)Q2{h)dhj 

as n 00. Since |/ a^{h; c)Qi{h)dh\ + \J a^{{0,h);c)Q2{h)dh\ + |cr2(0;c)| = 
0(1) as c -> 00, then by (l30l),(ini) and ([32]) . 

P[ lim limsupi;5f (c)/(ii's^^) = ) = 1. 
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Now we apply a classical Bernstein blocking technique for the proof of asymp- 
totic normality. Notations for the blocking technique of Bernstein are same with 
those of Lahiri. Let {Ai„} and {A2n} be two sequences satisfying the condition 
(A'6) and {A3„} = {Ai„} + {X2n}- Then the partition of the region ii„ is 
denoted by 

r„(/;e) = 7i(ei) X •••Jd(ed), e = (ei,--- ,ed)'e{l,2}^ 

where Ij{ej) = {IjX^n.lj^Sn + Ai„], if ej = 1 and Ij{ej) = {ljX3n + Ai„, (Zj + 
1)A3„,], if ej = 2. Note that with g(e) = [{1 < j < d : ej = 1}], 

for all I and e. Let eo = (1, • • • , 1)'. Then 

|r„(/;e)| =o(|r„(/;eo)|). 

Let Lin = {l '■ ^n{l', 0) C Rn} be the index set of all hypercubes r„(Z; 0) that 
are contained in R„, and let L2„ = {I : Tn{l; 0) D R„ 0, r„(?; 0) D R'^ ^ 0} 
be the index set of boundary hypercubes. With the notation above, S}^ can be 
separated into the sum of big blocks and small blocks and the sum of remaining 
variables. Here we consider only the case that station elements i and j are in 
the same block. If sums of pair whose elements are in different block, the joint 
probability of exceeding over the threshold would be zero as the sampling region 
is growing. Thus as n ^ oo, sums of pair would converge to and it could be 
negligible in consideration of our sum of processes. 

K 

Sk/<^k = Y^Wk{s'')r)k{s'')jaK 
fe=i 

|Li„| |Li„| 

q=l keJq 9=1 keHg feel-2n 

9=1 9=1 keL2n 

(big blocks + little blocks + leftover) 
where a% = Var( X;f=i ^if (sfe)6(sfe)) , •S'lx, = T,keJ,^kVk/(TK and 6"^^^ = 

J2keH,'^kVk/(TK- 

Two big blocks T{li; cq) and T{l2; eo) arc separated by the distance 
d{T{h; eo), r(Z2; eo)) > [(|Zi - Z2I - d)+X3n] + A2„. 
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By the strong mixing condition, 



leLi 



< C|Li„|a(A2„;A^). 



Therefore the asymptotic behavior can be shown with the independence of <S'^^^. 
Using Lemma A.l in |15) . we show that with probabihty one, 



9=1 

VariS'.^aK) = o{K' X'^^ s^^) , 



Now we have to show that 



J2ES[^;al-<rl=o{K'X-''sl^). 

9=1 

To prove above equation, we use Lemma 5.1 in ^3] and ([55)) - p5p . 



(33) 

(34) 
(35) 



(36) 



9=1 



< 



9=1 



< C J2C^M^iXltn^X-''' + \ognfami " h] ^ d)+X,,,] + X2n; X^) 

I /T^2\-2d 2 N 
+ 0{K A„ Si^) 

< C{d,Co)Ml{X\in^X--"' + \ognf{Xn/X^r.f''y^ 

A„/A3„ . 

a(A2„;At)+ ^ A:'^-ia(fcA3„ + A2„; A^J + o(i^2^-2S2^) 



fc=i 



o(A-A--.^^) 



Thus we show that the equation (p6|) holds and it is needed only to establish 
the Lindeberg condition. 



9=1 



0, as n —> oo. 
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Since we have 

\Li,.\ , 



/ I \ (2+<5)/2 . 

w) J\ 



<C\L,Jn-''y- 



< c 



2 + s \ ( 1 



g=l - I ^fceJ, ■wkml'^K\>(. 
\ (2+<5)/2 



2+<5 



keJg 

2+5 



2+5 



dP 



keJg 



dP 



'i^k^J„ ■<J'kVk/cK\>t 



keJg 



2+5 



dP 



— > as n — >■ oo, 
this imphes that the Lindeberg condition holds. 
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